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A class of discrete flavor-symmetry-based models predicts constrained neutrino mass matrix 
schemes that lead to speciflc neutrino mass sum-rules (MSR). We show how these theories may 
constrain the absolute scale of neutrino mass, leading in most of the cases to a lower bound on the 
neutrinoless double beta decay effective amplitude. 
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The discovery of oscillations [IHH] implies non- vanishing neutrino masses and mixing providing one of the most 
solid indications for physics beyond the Standard Model. The fact that neutrinos have very tiny masses, in contrast 
to charged leptons and quarks, and that two of the mixing angles are large, are among the deepest theoretical puzzles 
in particle physics. Since neutrinos carry no electric charge, they are expected on general grounds to be Majorana 
particles [9 , leading to the existence of lepton number violating processes [101 HI] • This intriguing possibility will be 
hopefully confirmed by the observation of neutrinoless double beta decay (Oz^/3/3) processes [12l[l3]. Indeed, upcoming 
QvjSp experiments are expected to improve the sensitivity by up to about one order of magnitude [TIHTT] . 

It seems unlikely that the observed pattern of lepton mixing angles is an accident: it probably indicates the existence 
of an underlying flavor symmetry of some sort, either an Abelian symmetry |18j or a non-Abelian one |19) . In the 
former case one typically obtains texture zeros for the mass matrices |20H22) but is unable to predict mixing angles. 
In contrast, non-Abelian flavor symmetries are potentially more powerful, allowing also in principle for mixing angle 
predictions. As an example, several realizations of non-Abelian discrete flavor symmetry schemes lead to an effective 
neutrino mass matrix which corresponds to a numerical (parameter-free) prediction for lepton mixing. A popular 
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example of such neutrino mass matrix M'^ is the tri-bimaximal (TBM) f32J type ^ , characterized by 



M"" = Mtbm 



^ X y y 

y X + z y ^ z 

\y y 



(1) 



z X + z 



which depends only on three complex parameters x, y and z. In the mass eigenstate basis the three complex parameters 
X, y and z would correspond to three neutrino mass parameters plus two Majorana CP phases [9l[70], as the Dirac 
phase disappears since = 0. Many such schemes are characterized by a specific (complex) relation among the 
parameters x, y and z f55H55] . leaving only two free complex parameters, further reducing the number of independent 
model parameters describing the lepton sector. In the mass basis these correspond to only two independent neutrino 
mass eigenvalues (the other follows from the existence of a neutrino mass sum rule), plus two Majorana CP phases 
(as mentioned, the Dirac phase is unphysical). 

In this paper we study the implications of these sum-rule schemes for the lower bound on the parameter |mee| 
characterizing the amplitude for neutrinoless double beta decay. We show that, given that two neutrino mass squared 
splittings are well-determined by neutrino oscillation data [H|S], we are left, approximately, with a one-parameter 
family of neutrinoless double beta decay theories in which the corresponding amplitude is mainly determined just 
by the overall absolute neutrino mass scale. The paper is organized as follows: in Section II we present the mass 
relations; in section III, we obtain the lower limit on \mee\ for all models considered here and briefly discuss their 
phenomenological implications, whereas in Section IV we present our conclusions. 



II. MASS RELATIONS 

In this section we focus on a general sub-class of mass matrices leading to a numerical (parameter-free) prediction 
for the lepton mixing matrix, consistent with current neutrino oscillation data [H [5] where the following types of mass 
relations hold: 



TTio mV ' 



B) 



C) X- 



(2) 
(3) 
(4) 



Here = m° denote neutrino mass eigenvalues, up to a Majorana phase factor, while x ^^id ^ are free parameters 
which specify the model, taken to be positive without loss of generality. For the sake of completeness, we also consider 
a fourth mass relation, 




As far as we can tell, this last relation has not yet been considered in the literature. In the following we show how 
this class of mass matrices arises in non-Abelian discrete flavor symmetry schemes. We first consider an effective 



^ There are also non-Abelian discrete flavor symmetry schemes leading to the bi- maximal lepton mixing pattern I23tf27l . as well as 
golden-ratio schemes |28H31| . 
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dimension-five operator description [7T!, and then discuss the cases where the neutrino mass matrix Ad^ arises from 
various seesaw mechanism reahzations, such as type-I or type-II seesaw [9l l72l - f77] or, from alternative low-scale seesaw 
schemes, for example, the inverse seesaw [75H5T| . 

Effective dimension-five operator description 

Consider first the dimension five operators (LLHH) , where the parentheses indicate all possible contractions among 
the irreducible representations of the underlying unspecified non-Abelian flavor symmetry group L and H belong to. 
Since we want a mass matrix with only two independent parameters, we assume that our effective Lagrangian contains 
only two terms, associated to two independent field contractions: 

£ = y^^e^L,L,),HH + ^^e\^{L.,L,\HH- (6) 

here M denotes the effective scale, a, h represent the two contractions, while e°- and are Clebsch-Gordan (CG) 
coefficients involving relevant field components Lij. After electroweak symmetry breaking, the effective neutrino mass 
matrix elements are given as linear combinations involving only the two parameters o and 6, 

Af^ = ae°+6e^, (7) 

where a = ya {Hf /M and h = {Hf /M. 

A number of non-Abelian discrete flavor symmetry realizations lead to the TBM structure for the effective neutrino 
mass matrix in Eq. ([ij with some suitable relation among the x, y and z coefficients. It is clear that in this case 
the corresponding mass eigenvalues will always be expressed as linear a combination of a and h and hence will be 
related to each other through a relation of type {A) in Eq. ([2]). 

Type-I seesaw mechanism 

For the type-I seesaw mechanism there are two simple ways to get a neutrino mass matrix similar to M^bm 
depending only upon two free complex parameters. In the first case, the Dirac neutrino mass matrix m^j has the 
structure given in eq. ([!]) while the right-handed neutrino mass matrix Mr is proportional to a numerical fJ- — t 
invariant matrix satisfying the relation (2, 2) + (2, 3) = (1, 1) + (1, 2) among its elements, like for instance in Eq. (|8|: 

/lOO\ / 2 I I \ 

Mfloc 1 , -11-1, 15-2,... (8) 
\001/ \-l-ll/ \l-2 5/ 

In what follows we call such matrix generically as "TBM-type" . It is not difficult to verify that the light neutrino mass 
matrix arising from the type-I seesaw formula [82 AI'^ = —raj^M^^m^ has mass eigenvalues of type m}^ cx (aia-|-/3i&)^, 
yielding mass relations of type (C). For instance, in Ref. [83] it has been found that m\ cx (a + 6)^, oc a? and 
TOg (X {a — b)^, from which the relation ^/m^ + -yAi^f = ^^m!^ holds. 

The second possibility arises when M^, ^ Mtbm as in Eq. ([T]), while the Dirac neutrino mass matrix is a numerical 
"TBM-type" matrix, as in eq. (|8]). In this case it is simple to show that the eigenvalues of are of the form 
m\ cx l/{aia -t- f^ih), where and Pi are numerical coefficients, giving a mass relation of type (B). For instance, in 
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the model of Ref. [31] the authors found rrii cx l/(a + &), m2 oc 1/a and cx 1/(6 — a) from which the relation 
l/rrii — = 2/m2 is satisfied. 

Other seesaw mechanisms 

Similar conclusions can be obtained for different seesaw mechanisms, such as type-II. From the point of view of our 
classification, type-II seesaw is equivalent to the dimension five operator case (A). 

We now move to the inverse seesaw mechanism |78| [79] , which arises when introducing a fermion singlet S with 
opposite lepton number with respect to the right-handed neutrinos, so that the effective light neutrino mass matrix 
is nil, = 'tzd jj/i jjm|^. Assuming mn and /z to be proportional to the identity matrix and M ~ Mtbm, it is 
straightforward to show that we can obtain the mass sum-rule of type (D). 

A novel seesaw mechanism arises from left-right symmetry [ST] or the full SO(IO) [M] in the presence of gauge 
singlet fermions, and has been called the linear seesaw. In such scheme the effective light neutrino mass matrix is 
given in terms of three independent sub-matrices and scales linearly with respect to the usual Dirac neutrino Yukawa 
couplings, hence the name. In order to have a mass relation, we need two sub-matrices of numerical "TBM-type" 
like in eq.([8]) and the third one similar to Mtbm (see also [HS]): otherwise additional free parameters are introduced, 
beyond our assumed two. One can show that all four cases can be realized, depending on which matrix has the form 
Mtbm- 

III. LOWER BOUND FOR NEUTRINOLESS DOUBLE-^ DECAY 

Let us first consider the amplitude for neutrinoless double-/? decay within a flavor-generic model. One can plot the 
effective neutrino mass parameter jmeel determining the Oi^/3/3 decay amplitude, as a function of the lightest neutrino 
mass. As is well-known, by varying the neutrino oscillation parameters Arn^^^, Am^^;, 9i2, 0i3, 023 in their allowed 
ranges [H[5] one obtains two types of relatively broad bands in the ( | rriee |, wfj^j^j) plane corresponding to normal and 
inverse hierarchy spectra, as shown in Fig. [l] 




FIG. 1: Comparison on the allowed range of (|mee|) as a function of the lightest neutrino mass. For the TBM mixing pattern 
(red and green bands for NH and IH respectively) and for the full allowed Scr C.L. ranges of oscillation parameters from [31 [5] 
(gray and blue bands for NH and IH respectively). 

In this "generic" case there is a lower bound on the neutrinoless double-/? decay effective mass parameter |mee| only 
in the case of inverse mass hierarchy: due to the possibility of destructive interference among the light neutrinos from 
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the effect of having opposite CP signs or due to the effect of Majorana phases, no lower bound can be estabhshed for 
the case of normal hierarchy [85H87] . 

Let us now turn to the case where MSR relations like (A),(B),(C) and (D) hold. As discussed above these can be 
obtained in flavor models where the neutrino mass matrix only depends on two independent free parameters, so that 
the resulting mixing angles are fixed, like for example for the tri-bimaximal or bimaximal mixing patterns. 

For definiteness here we focus on the case where the rotation in the neutrino sector is of tri-bimaximal form. 
Corrections from higher dimensional operators and/or from the charged lepton sector can yield ^13 7^ 0, as suggested 
after the T2K 12] and Double-Chooz [3 first results [4 . 

Hence we retain the TBM approximation as a useful starting point to obtain our MSR relations. However, when 
evaluating a lower bound on the effective neutrino mass parameter |mee| determining the neutrinoless double- /3 decay 
amplitude, we include explicitly the effects of non- vanishing ^13. We do this by taking the values at 3 cr determined 
in Ref. Such a lower bound can be obtained from the following procedure. 

We first consider that the neutrino masses are complex parameters, where the two Majorana phases are encoded 
in m2 and TO3, i.e. 



Ad 



nVie 



(9) 
(10) 
(11) 



As shown in Fig. [2] the neutrino mass sum-rule can then be interpreted geometrically as a triangle in the complex 
plane, whose area provides a measure of Majorana CP violation ^. Each of the above equations (|9|), (10) and (11) 




FIG. 2: Neutrino mass sum-rule triangle in the complex plane. 

can be split into two independent equations for the real and imaginary parts. 

For simplicity let us start from the idealized case where the neutrino oscillation parameters Arn^j^, Am^^;, 6'i2, 
^13; ^23 are perfectly well-measured quantities. In this case one can extract the two Majorana phases a and /3 as 
functions of the base of the triangle, which is determined by in case of normal hierarchy (NH) or by mg in case 
of inverted hierarchy (IH), as well as the parameters x and f labeling the particular model under consideration. 



^ As the area shrinks to zero one obtains the CP-conserving hmits corresponding to the four independent choices of CP sign |85l I86| . 
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These relations obtained can then be inserted into the general expression of |rneeh 

I'TT-eel = |ci2Ci3TOl + S^2Ci3e*"m2 + S^ge^'^TOa | . (12) 

For each (x, model this effective mass parameter depends on a single parameter, namely the length of the triangle 
base, which gives a measure of the absolute scale of neutrino mass. 
For instance, for case (A) this procedure gives: 

cos a m2 + X cos /? TTig = TO? , ^sina + x sin/3 mg = , (13) 

so that the Majorana CP phases are determined as: 



2 ^,2 i A^2 



cos a 



cos/3 — 



2Wl^-\/ Ato^^; + TO 

{ml + AmltJ - e {r 
2toiXv/Ato2j„ +TO 



2 



(14) 



The lower bound for the lightest neutrino mass can be obtained from our MSR, using the triangle inequality in the 
complex plane as suggested by Rodejohann and Barry in : 88] . see Fig. ([2]) for a schematic view. 

We must first select the biggest side of the triangle; calling them xi, X2 and 0:3, then the triangle inequality 
\xi\ < \xj\ + \xk\ must be fulfilled, where \xi\ = Max(|a:;i|, \x2\, [xaD and i ^ j ^ . In case (A) and assuming NH for 
the neutrino mass spectrum, we always have (xto2,^TO3) > to? and the largest side of the triangle can be either XTO2 
or ^TO§; so we must consider separately these two cases. After rewriting two masses in terms of the two squared mass 
differences, we can obtain a lower limit for the lightest neutrino mass from the triangle inequality \xi\ < \xj\ + \xk\- 
For the other cases we follow the same procedure. The lower bound on the lightest neutrino mass obtained in this 



way is then used to estimate the lower bound for |TOee| from the general expression in eq. ( 12 ). Notice that, although 
we focus here on TBM schemes, some of the MSR considered in our analysis can also be derived using bimaximal 
mixing as a starting point: for instance, in Refs. [26l[27] a relation of type (A) with (x,<^) — (lj2) has been derived 
and the phenomenological consequences studied. 



Classification 

The results obtained from the procedure discussed above are summarized in Tabjlj where we report the lower limits 
of |TOee| corresponding to different integer choices of {x,0 between 1 and 3 and for each of the four MSR considered 

. The cases already discussed in 
the literature are indicated by giving the corresponding reference. The entries denoted with the symbol (★) represent 
situations that do not satisfy the inequality for any value of the lightest neutrino mass rriught- Cases marked by a (-) 
correspond to models which, as far as we can tell, have not been considered. Some comments are in order. First let 
us consider the effect of a possible non-zero effect of 613 as indicated by recent experiments [3] as well as global 
neutrino oscillation fits [3H1]- In Figlslwe show the prediction for \mee\ as function of mught obtained from the MSR 



in eqs. ([2])-(l5]), for both normal and inverted hierarchies (compare also with Figs. 4]l7 



^ Notice that there are three inequahties of the type \xi\ < \xj \ + \xk \ obtained by permuting the three indices i, j and k, but only one of 
these constrains the hghtest neutrino mass. 
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A-NH 


A-IH 


Rcf. 


B NH 


B - IH 


Rcf. 


C- NH 


C -IH 


Rcf. 


D -NH 


D- IH 


1 1 


0.010 


0.044 


,33r.49j 


0.008 


0.036 


50rl5'51 


0.006 


0.029 




.005 


.008 


1,2 




0.046 




0.008 


0.027 




■A- 


0.014 




.004 


.026 


1 3 




0.011 




0.030 


0.005 






0.014 




.018 


.025 


2 1 


0.006 


•A- 


[41, 53 


0.006 


0.007 


[49, 54 68 


0.000 




IS3I 




.007 


2,2 


0.019 


0.026 




0.023 


0.008 




0.017 


7k- 




.003 


.015 


2,3 


★ 


0.046 




0.007 


0.008 






0.031 




.005 


.026 


3,1 


0.004 


★ 




0.004 


0.008 




★ 






★ 


★ 


3,2 


0.011 






0.004 


0.021 




0.000 








.007 


3,3 


0.023 


0.061 




0.029 


0.031 




0.011 


0.019 




.018 


.016 



TABLE I: Minimal values for the effective Qv^P decay mass parameter |mee|, in eV, see text for details. 



3^7712 + 3-^7713 = (right panel) and 2^TO2 + ^/^s ~ \/™T (left panel). For the red bands we assumed the TBM 

values of the oscillation parameters (implying 6*13 = 0) while in the yellow bands corresponds to the same MSR, but 
now varying the values of ^13, ^23 and 9i2 within their 3cr C.L. interval. By looking at the left panel in Fig. [3] one 
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FIG. 3: {|mee|) as a function of the lightest neutrino mass corresponding to the mass sum-rule 2y/m2 + y/rrMi ~ ^rn\ '83' (left) 
and 3.^777,2 + Jt^fmi — (right). The red bands correspond to the TBM mixing pattern, while the yellow bands correspond 

to the same MSR, but now varying the values of 6113, S23 and 612 to 3(7 C.L. range. 



sees that, indeed, the 07//3/3 lower bound is sensitive to the value of ^13. 

One also finds that, as expected on general grounds, all inverse hierarchy schemes corresponding to various choices 
of (x,C) within sum-rules A-D have a lower bound for the parameter |?77ee|- However, the numerical value obtained 
depends on the MSR scheme, signaling that not all values within the corresponding band in Fig. [l] are covered. 

On the other hand, even though normal hierarchy models do not lead to a lower bound on the OvPfi amplitude 
due to the possibility of destructive interference amongst the light neutrinos, one finds that the possibility of full 
cancellation is precluded for all schemes in the table, except for the (2,1) case considered in Ref. [SJ and the (3,2) 
scheme, both of which correspond to MSR of type (C). All other NH MSR schemes considered here imply a minimum 
value for the Oiy/S/S decay amplitude One finds that the most favorable cases are given by: 



* Of course some of the bounds are phenomenologically less interesting since they fall outside reahstic sensitivities of coming experiments. 



8 



{x,0 = 


(3,3) 


for the 


case 


(A) NH, 


(x,0 = 


(3,3) 


for the 


case 


(A) IH, 


(x,0 = 


(1,3) 


for the 


case 


(B) NH, 


ix,0 = 


(3,3) 


for the 


case 


(B) IH, 


(x,0 = 


(3,3) 


for the 


case 


(C) NH, 


(x,0 = 


(2,3) 


for the 


case 


(C) IH, 


(x,C) = 


(3,3) 


for the 


case 


(D) NH, 


(x,C) = 


(1,2) 


for the 


case 


(D) IH. 



In particular, the maximal value we have found for the lower bound on |mee| is \mee\ = 0.061 eV, obtained in 
correspondence with the set of values (x,C) — (3,3) for the case (A) in IH. Such a value for \mee\ lies within the 
sensitivity of upcoming experiments; hence it would be interesting, from the model building point of view, to find 
from first principles a flavor-symmetry-based model predicting such a mass relation; we will return to this problem 
elsewhere. 

The same phenomenologically interesting cases are now studied more in detail, showing the behavior of \mee\ as 
function of the lightest neutrino mass in Figs. |4|5|6| and [7] In all plots, the two bands are the most generic ones 
compatible with both normal and inverted hierarchies, derived considering the 3cr allowed ranges on the neutrino 
oscillation parameters as obtained in Ref.p^ and consistent with latest T2K and Double-Chooz experiments, see Fig. 

m 

In Fig.|4]we give the allowed (|TOee|) values as a function of the lightest neutrino mass. The figures correspond to case 
(A). In the left panel the yellow band corresponds to the model which predicts the mass sum rule 3m2 -I- 377i3 — nii 
in case of NH. On the right the red band corresponds to the same sum rule in the case of IH. Other MSR Oi'/3f3 
amplitude lower bounds are illustrated in subsequent figures. 
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FIG. 4: Effective OvfiP mass parameter {|mee|) as a function of the lightest neutrino mass for MSR type-A schemes. On the 
left panel the yellow band corresponds to the model which predicts the sum mass rule 3m2 -I- Sms — mi for the case of NH. 
On the right the red band corresponds to the same sum rule in the case of IH. 



IV. CONCLUSIONS 

In this paper we have analyzed the implications for the lower bound on the effective 0i>/3f5 neutrino mass parameter 
|mee| arising from possible mass sum-rules obtained in the context of fiavor models. Mass sum rules are classified in 
four different categories, some have already been considered in the literature. For each case, we have first extracted 
the allowed numerical values of |mee|, for both mass orderings of the neutrino mass eigenstates and we have then given 
the behavior of |mee| as a function of the lightest neutrino mass. Although our MSR schemes were obtained within 
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m„(eV) m„(eV) 

FIG. 5: (|mee|) as a function of the lightest mass. The figures correspond to the case (B). On the left panel the yellow band 
corresponds to the model which predicts the sum mass rule l/m2 + S/ms = 1/mi for the NH case. On the right, the red band 
corresponds to the prediction of the sum rule 3/m2 + S/ma — 1/mi in the case of IH. 




m„(eV) m„(eV) 

FIG. 6: {|mee|) as a function of the lightest neutrino mass. The figures correspond to the case (C). On the left the yellow band 
corresponds to the model which predicts the sum mass rule 3^/m2 + S^/m-j = ^/rnl in case of NH. On the right, the red band 
corresponds to the prediction of the sum mass rule 2^m2 + S^ma — ^mi in the case of IH. 




m„(eV) m„(eV) 

FIG. 7: (|mee|) as a function of the lightest neutrino mass, case (D). On the left the yellow band corresponds to the MSR 
scheme which predicts 3/^m2 + S/^ms = l/^/mi for NH. On the right panel, the red band corresponds to the prediction of 
the mass sum rule l/^m2 + 2/^1113 = l/y^mT in the case of IH. 



the TBM anzatz, we have computed the possible values of |mee| considering all the neutrino parameters (including 
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a non- vanishing ^13) within their 3a allowed ranges. In most MSR schemes one finds a lower bound for the 0^/3(3 
amplitude, even for NH spectra. We find that the most favorable case (large lower bound) corresponds to a sum-rule 
of type (A) obtained in correspondence of the set of values {x,£.) = (3,3), |mee| = 0.061 eV. Such a mass relation has 
not been considered so far, and the searching of a flavor model able to predict it at leading order is now in progress. 
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